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Abstract. We have investigated the heaviest one-neutron halo candidate ^^C 
nucleus. Few-body model calculations of cross section angular distributions for 
the ^^C(p, d)^*C reaction, together with the test calculations carried out for 
the C{p^dY^C reaction, at a low incident energy are presented for different 
possible halo-neutron configurations. We show that there is a clear distinction 
between in particular In = and In = 2 halo transfers. The sensitivity of the 
cross sections to the assumed ^^C single neutron separation energy is discussed. 

1 Introduction 

Dripline physics is at the focus of present-day nuclear science. Experimentally, 
thanks to existing and emerging radioactive-ion-beam facilities, we are on the 
verge of invading the territory of extreme N/Z ratios in an unprecedented way. 
Theoretically, nuclear exotica represent a formidable challange for the nuclear 
many-body theories and their power to predict nuclear properties far from 
stability. In order to give a unified description of all nuclei in the periodic table, 
theoretical systematic studies in a single model are desired for both exotic nuclei 
near the particle-drip line and ordinary nuclei near the stable line. Certainly, 
this is a heavy task for future research. 

A series of experiments have shown that there exist neutron halos in ex- 
otic light nuclei near the neutron-drip line. These studies change our views on 
nuclear structure and lead to a new impact to the traditional nuclear physics 
along the isospin freedom. Going to the limits of the nuclear binding is also 
important for an improvement of our description of normal nuclei from the 
neighborhood of the beta stability valley. As our knowledge on nuclei far from 
stability deepens, it becomes more apparent that nuclear halos are a general 
feature of loosely bound nuclei. As the binding energy becomes smaller in the 
vicinity of the drip lines, the valence nucleon(s) tunnel out of the central poten- 
tial, and enhance the diffuseness of the nuclear surface. Eventually, this leads 
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to a delocalization of the valence nucleon(s) which can be pictured as a halo 
surrounding the core of the nucleus. The appearance of this halo is determined 
by the height of the potential barrier, which itself depends on the binding 
energy, the angular momentum and, for protons, the Coulomb potential. Ex- 
perimental and theoretical studies of exotic light nuclei with a normal, localized 
nuclear core and a dilute few-neutron halo or skin are now well advanced and 
are becoming increasingly sophisticated. Such systems provide a stringent test 
of nuclear structure models developed for stable nuclei as they involve new 
structures and surface phenomena. 

The question about the halo nature of ^^C and its underlying structure 
is one of the interesting current questions in dripline physics. The neutron 
rich ^^C nucleus, the last particle-stable odd-neutron isotope of carbon, has 
drawn much attention recently as candidate for having a one-neutron halo 
structure This can be easily anticipated due to its very small neutron 

binding energy (of the order of a few hundred keV). However, halo formation 
not only depends upon the small separation energy of the valence neutron. 
The spin-parity and the configuration of the single particle state of the valence 
neutron also play important roles in the formation of the halo. The case in 
which ^^C has a ground state = 1/2+ with a dominant ^^C (0"'")(8>lsi/2 
single particle configuration favours the halo formation. On the other hand, 
the case of a J'^ = 5/2+, expected from standard shell model ordering, prevents 
the halo from being formed due to the large centrifugal barrier associated with 
a 0^5/2 orbit. The — 5/2+ may also result when the dominant single particle 
configuration in the ground state is ^^C (2+) ® lsi/2- In this case, the binding 
energy of the valence neutron is effectively increased by the excitation energy 
of 1.62 MeV of the excited (2+) state of the core, thereby also hindering 
the halo formation. A — 3/2+ is also possible with this configuration for 
the ground state, as considered in the recent model calculations 

In order to probe the structure of ^^C, only a few experiments on the 
dissociation of this nucleus have been carried out so far Very 
recently the relative energy spectrum of ^^C has been measured in its Coulomb 
dissociation at 67 MeV/ nucleon incident energy on Pb target in a kinematically 
complete experiment at Riken 1^ . The Riken data, when compared with semi- 
classical calculations H], support a dominant lsi/2 configuration of ^^C halo 
wavefunction. The MSU data on the other hand, suggest an s— wave neutron 
around the 2+ excited state of ^^C. Recent theoretical calculations on the 
Coulomb breakup of ^^C ||] indicate the possibility of an s— wave neutron 
coupled to the 2+ excited state of ^*C, or a d— wave neutron, or a coherent 
superposition of these two configurations as far as the recent Ganil data Q 
are concerned. The total one-neutron removal cross sections in the Coulomb 
dissociation of ^^C on heavy targets have been calculated in |^. These cross 
sections, when compared with experimental data one-neutron removal cross 
sections measured at Ganil [Q and MSU ||^, show that the configuration of a 
lsi/2 neutron coupled to the 0+ ^^C core has non-significant contribution to 
the ground state structure of ^^C ^. In view of all these discrepancies, it has 
not been possible to draw any definite conclusion on the structure of ^^C. 
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The difRculties in drawing conclusion about the structure of ^^C might arise 
also due to lack of precise knowledge of the one-neutron seperation energy 5„ 
in ^®C. This quantity is very crucial as far as the property of the neutron halo 
is concerned: the asymptotic behaviour of the radial wavefunction is governed 
by the separation energy and the extension of the halo is therefore directly 
linked this quantitiy. The value used in most of the calculations to date is 
240 keV. This value is a world wide average of several experiments each with 
significant uncertainties pO{ . The recent measurement of the relative energy 
spectrum of ^^C in its breakup on Pb at Riken has led to an indirect determi- 
nation of the ^^C one-neutron separation energy 1^. A semi-classical analysis 
of this new, kinematically complete, measurements of the Coulomb dissociation 
of ^^C showed the data to be consistent with 5„ « 530 keV. In the present cal- 
culations, we will consider the consistency of this proposal with the interaction 
cross sections. 

In this brief report, we have re- investigated the structure of ^^C through 
the study of pickup reactions involving ^^C nucleus at low energies, as recently 
the use of low energy single nucleon transfer reactions for structure studies 
of exotic nuclei have attracted attention Because of the simplicity of 

the theoretical interpretation of such reactions, they are thought to provide an 
important source of information about the structure of halo nuclei. 

Single nucleon transfer reactions, such as the {p, d) and {d,p) reactions, have 
been a reliable tool in nuclear spectroscopic studies of stable nuclei. Besides 
investigating the energy level structure of such nuclei, the spins and parity 
(J^) of levels can be determined from the shape of the differential cross section 
angular distribution. The shape of an angular distribution of an unknown level 
is either compared to those of known levels (ideally in the same nucleus), or to 
Distorted Wave Born Approximation (DWBA) calculations. By this technique, 
at least the transferred orbital angular momentum {£„) may be determined, 
and for J'^ levels with relatively simple structure, one can usually then deduce 
J'^ from shell model arguments. 

Since radioactive beams of sufficient intensity are now starting to become 
available, low energy single nucleon transfer reactions are again being proposed 
as a spectroscopic tool, this time in inverse kinematics, since /3— unstable nuclei 
are produced as secondary beams, with a proton or deuteron target for the 
radioactive beam to investigate the structure of nuclei far off stability. Since 
single nucleon transfer reactions are binary processes the interpretation of the 
results is comparatively simple and decisive answers on nuclear structure and 
also insight into the reaction dynamics of exotic nuclei can be expected. The 
discussion in this report is limited to the {p, d) reaction. A more general review 
of light-ion transfer reactions in inverse kinematics may be found in Ref. [l^ . 

It would thus be of considerable interest to measure the spin and parity 
of ^^C via the (p, d) reaction. The shape of the angular distribution of the 
^'^C{p,dY^C reaction, measured in inverse kinematics, might be used to dis- 
tiguish between £„ = and = 2 transfer, if the last neutron is in a pure 
single particle state. Unfortunately, to the best of our knowledge, the beam 
intensity of ^^C produced by present facilities put this out of reach for the 



4 



moments. Nevertheless, the resuhs of the present work would be beneficial for 
the experimentalists in the near future. 

A measurement of the angular distribution, hovewer, for the ^'^C{p,d)^^C 
reaction is more feasible. However, like ^^C, the ground state spin and parity 
for ^^C is not known experimentally. Again, the energies of the 5/2+, 3/2+, and 
1/2+ levels given by the shell model are too close to constitute a prediction, 
although a comparison of the experimental and theoretical beta-decay lifetimes 
favours 3/2+or 5/2+. Knowledge of the ^^C ground-state spin would help de- 
termine the shell model interaction in this {A, Z) region, which would lead to 
better confidence in the predictions for ^^C. Hence, we start first with the anal- 
ysis of the ^'^C{p, dy^C reaction before investigating the cross section angular 
distributions of the ^^C{p, d^^C reaction for different possible halo- neutron 
configurations. 



2 Application 

The most dramatic feature of halo nuclei is their very loosely bound few-body 
character, with a strong spatial localization of the core nucleons and a delo- 
calization of the halo particle(s). More recent theoretical analyses have shown 
that an explicit treatment of this correlated few-body nature, or granularity, is 
important quantitatively for calculations of reaction cross sections . Such a 
few-body description leads to larger/smaller calculated reaction cross sections 
depending upon the scattering angle than those obtained from the standard 
calculations which do not consider the loosely bound structure. Such descrip- 
tion has significant implications for the deduced size and ground state structure 
of the halo. It follows that interaction cross section analyses are model depen- 
dent and require the use of theoretical relative motion wavefunctions for these 
few-body structures. 

We discuss the exit channel of the ^^C(p, d^^C and ^^C{p, d^^C reactions 
within three-body (halo neutron (n) + p + core nucleus) models. We calculate 
the transfer amplitude using the prior form of the (p, d) matrix element, thus 
the transition interaction is the n—p interaction and we need a full (three-body) 
description of the n-\-p-\-core system in the final state. For the description of this 



final state we have used the Adiabatic (AD) model 1 15 and the Quasi-adiabatic 
(QAD) approach To clarify the importance of the breakup corrections to 
the transfer cross sections, we also perform the standard DWBA which do not 
account for the effects arising from the breakup of the projectile in the field of 
the nucleus. 

In the AD treatment it is assumed that the excitation of the projectile is 
to states in the low energy continuum; its treatment of possible high energy 
breakup contributions is therefore naturally suspect. An improved treatment 
of these higher energy breakup configurations is provided by the QAD method 
calculations which takes approximate account of modifications to the centre of 
mass energy of the n — p pair in breakup configurations through the use of a 
mean breakup relative energy for the continuum states. It thus breaks the de- 
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generacy with the elastic channel. The inclusion of these higher energy breakup 
components leads, in general, to an enhanced transition amplitude and a better 



description of experimental data |16 . These quantum mechanical three-body 
models have been well described in Ref. by the work of Yilmaz and Goniil 
through a successful application of the models in analyzing a halo transfer 
^^Be{p,d)^'^Be reaction. The details of the resulting partial wave expansions 
and solution of the equations can also be found in Sec. Ill of Ref. [ p^ . 

The calculated two-body (DWBA) wavefunction, using deuteron folded po- 
tential, and three-body (AD, QAD) wavefunctions using the same inputs, in 
partial wave form, are employed in a modified version of the program TWOFNR 
for the evaluation of the reaction observables, performing the zero-range 
approximation. The radial integrals are carried out from to 40 fm in steps 
of 0.1 fm. The maximum number of partial waves used is 30 for both entrance 
and exit channels. As the long tail of the halo wavefunction makes internal 
contributions less important, all calculations presented in this paper are done 
without non-locality corrections. 



2.1 The ^^C{p, df^C Reaction at 30 MeV 

As an additional test, we first deal with the ground state of ^^C. The separation 
energy of the last neutron in this nucleus is again small, only Sn ~ 730 keV [ p^ . 
Data concerning the structure of the ^^C nucleus is quite poor and can easily be 
summarized in a few sentences [|l9| . Observation of beta-delayed neutron emis- 
sion has been reported An excited state of ^^C is observed at 295 ±20 keV . 
Three closely spaced low- lying states are expected again like for the ^^C case, 
i.e. J'^ = 1/2+, 3/2+, 5/2+, and it is not clear which one of them is the ground 
state. Nevertheless, based on the systematics described in j2^], where a modi- 
fied Millener-Kurath interaction is employed within the shell model formalism, 
it is unlikely that ^^C has a ground state J"^ = 5/2+. Morever, analysis of the 
gamma-ray spectra following the beta-decay of ^^C does not support 5/2+ as 
the ground state spin either Unfortunately, as is clearly stated in [^ , 
the present shell model calculations cannot give a definitive prediction for the 
ground state spin of ^^C, and 1/2+ and 3/2+ both remain viable possibilities. 
On the other hand, both shell model ^ and relativistic mean field calculations 
[ p2| argue for a J'^ = 3/2+ non-halo ground state structure for this case. In 
a recent paper |p3| , neutron capture reaction rates on neutron- rich C— , N—^ 
and O— isotopes have been calculated in the framework of a hybrid compound 
and direct capture model. The results have been compared with those of pre- 
vious calculations as well as experimental results. In that work, Herndl et al 
have obtained the neutron spectroscopic factors for the possible ground states 
of ""^^C We make use of these spectroscopic factors in the present calculations 
discussed in this section for more realistic absolute values of the cross sections. 

To investigate the difference in shapes of the angular distributions with the 
various possible ground state spin assignments, we have performed calculations 
with the zero-range approximation. The binding potentials in all cases are taken 
to be of Woods-Saxon type with the radius and diffuseness parameters as 1.25 
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fm and 0.65 fm respectively. Their depths have been calculated to reproduce 
the corresponding binding energies. The entrance and exit channel potentials, 
including spin-orbit interactions, are obtained from the global parameterization 
of Bechetti and Greenless . 

We have found that the shapes of the angular distribution are not sensi- 
tive to the radius of the neutron binding potential. Test calculations involving 
different sets of potential parameters have shown that the exact location and 
sharpness of extrema in the cross sections depend on the choice of optical model 
potential but there is always a clear distinction between £„— transfers. In par- 
ticular, the £n = cross section always has a maximum at 6c.m. — 0° and falls 
off rapidly with angle. From the calculations shown in Figs. ?? and ??, it is seen 
that at 30 MeV/nucleon one would require to measure cross sections out to 
at least 9c.m. — 15° to distinguish clearly between £n — and = 2 transfer. 
Despite the uncertainties on optical potentials which in fact are important for 
a precise description of transfer reaction cross sections, the results can be ex- 
pected to describe realistically the dynamical features of such reactions. Thus, 
such a measurement, besides being of interest in its own right, would constitute 
a first step towards a future ^^C experiment. 

In addition. Fig. ?? makes clear that the inclusion of breakup effects in halo 
transfer reactions produces changes in the shape of the angular distributions 
and these effects increase the absolute values of the theoretical cross sections 
in forward direction and thus lead to the smaller values of the spectroscopic 
factors extracted from the experimental data. It is expected that such effects 
will be larger when the incident energy increases and the mass of the target 
decreases. We thus will only consider the more exact QAD model calculation 
results in the next section in treating the deuteron channel of the ^^C{p, d)^^C 
reaction. 

2.2 The ^^C{p, dy^C Reaction At 30 MeV 

Since the spin/parity for the ground state of ^^C is still experimentally un- 
known, we have considered in our calculations the following configurations for 
the valence neutron in ^^C: 

(a) a lsi/2 state bound to a 0+ ^^C core by 0.24 MeV, 

(b) a lsi/2 state bound to a 2+ ^^C core by 1.86 MeV, ^ 

(c) a 0^5/2 state bound to a 0+ ^^C core by 0.24 MeV, and 

(d) a lsi/2 state bound to a 0+ ^^C core by 0.53 MeV. 

Since the separation energy Sn is not known with satisfactory accuracy it 
is treated as an input parameter in our calculations. Therefore we have also 
performed calculations (d) for a larger value, Sn ~ 530 keV, to investigate the 
sensitivity of the results to the separation energy. 
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The parameterization procedure of the ^^C core and the single neutron 
interaction within the QAD model is again taken similar to the case of the pre- 
viuos calculations. We initially, see Fig. ??, assume a Woods-Saxon n + C 
binding interaction with radius parameter rp = 1.25 /to, diffuseness oq = 0.65 
/to, and neutron separation energies 0.24 MeV and 1.86 AleV for the core 
ground and excited state configurations, respectively. However, in order to 
consider the cross section sensitivity to the assumed neutron binding poten- 
tial geometry (ro, ao), we perform our calculations with different sets of poten- 
tial parameters. We see that the cross sections at large angles increase with 
increasing ao and ro parameters while they seem similar at forward angles. 

As seen from Fig. ??, reducing Sn does not considerably change the shape 
and magnitude of the cross section angular distributions for the reaction con- 
sidered, unlike the case (b) that leads to the larger cross sections with a different 
absolute value at forward angles when compared to the case of (a). 

In Fig. ??, the energy dependence of the cross section for a possible neutron- 
halo configuration of the reaction under study is shown. This observation might 
be of interest for the design of experimental setups. The small separation en- 
ergy and large spatial extension of the neutron-halo wavefunction have lead to 
an energy dependence of the cross sections which is very different from these 
known for transfer reactions on well bound systems. The rapid decrease of 
the calculated cross sections with increasing incedent energy for the weakly 
bound nucleus shows that the dynamics of transfer reactions is favorable for 
systems far off stability and experiments can very likely take advantange of 
this behaviour. Only transfer reactions with rather large cross sections can be 
investigated due to the low intensity of secondary beams and the use of in- 
verse kinematics in reactions involving halo nuclei. It is clear from Fig. ?? that 
low energy transfer reactions indeed might be useful for investigations of such 
short-lived isotopes. 

3 Summary and Conclusion 

Recently experimental evidence for a new one-neutron halo candidate, namely 
^^C, has become available. This has sparked off considerable interest as the 
only established one- neutron halo nucleus so far has been ^^Be where relative 
s— motion dominates the ground state. The ^^Be nucleus has been extensively 
studied during the last few years and has provided a testing ground for single 
neutron halo theories [12, and the references therein]. 

The structure of the ground state and its last neutron separation energy 
Sn are however still very uncertain. A naive shell model suggests a nodeless 
orbit for the least bound neutron. More detailed calculations |^ however 
predict a 1/2+ ^^C ground state due to a lowering of the lsi/2 orbital. A 3/2+ 
or 5/2+ ground state remains a theoretical possibility ||], but would involve 
(2+(X>lsi/2) and/or (2+ ® 0^5/2) excited ^^C core components, and (2+(g)lsi/2) 
and/or (0+ 00^5/2) components, respectively. 

To investigate this uncertainty, we have taken into account the idea that 
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in combination with independent experimental data, such as from Coulomb 
dissociation measurements and from momentum distributions following ^^C 
breakup reactions, cross section angular distribution measurements for halo 
transfer reactions involving ^^C nucleus may at this stage provide a useful con- 
straint on the ground state of such nuclei. Along this line, we have studied 
the ^^C{p,dy^C neutron-halo transfer reaction, together with the test cal- 
culations on the ^'^C{p,dY^C reaction, performing the calculations for both 
lsi/2 and 0(^5/2 neutron configurations with separation energy Sn = 0.24 MeV 
about a ^^C{0^; g.s) core, and with separation 1.86 MeV about 

a ^®C(2+; 1.62 MeV) excited core. The comparison between the calculation 
results, which have been obtained using approximate two- and three-body 
quantum mechanical models for different possible single neutron configura- 
tions, gives some support to the suggestion that at 30 MeV/nucleon incident 
energy one would require to measure cross sections out to at least Oc.m. = 15° 
to distinguish clearly between = and ^„ = 2 transfer for the reaction of 
interest. 

The (p, d) reaction is a potentially useful spectroscopic tool in inverse kine- 
matics with radioactive beams. Its advantages include its simplicity (reliability 
of model calculations) and relatively large cross sections for single particle 
states. A promising aspect for the future is that the instrumental energy and 
angular resolution requirements will be more relaxed. The high intensities of 
radiaoactive beams will not only allow transfer reactions to be performed with 
exotic isotopes, but will also permit coUimation of beams close to the stability 
line, thus giving better emittance quality. 
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FIGURE CAPTIONS 

Fig. ?? Calculated differential cross section angular distributions obtained 
by different theoretical models for the reaction ^^C(p, d)^^C at 30 MeV for 
a lsi/2 , 0^3/2 and 0^5/2 transfer. The spectroscopic factors are set to 0.644, 
0.035 and 0.701, respectively. 

Fig. ?? Calculated angular distributions using the QAD model for the tran- 
sitions as in Fig. ?? 

Fig. ?? Calculated differential cross section angular distributions obtained 
by the QAD model for the ^^C{p, dY^ reaction at 30 MeV for possible neutron- 
halo configurations. The spectroscopic factor is 1. 

Fig. ?? As for Fig. ??, but the QAD calculations have been carried out for 
different bombarding energies in case of a (0+ (g) lsi/2) neutron-halo configura- 
tion. 
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